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Charged-Higgs effects in a new B DtVt differential decay distribution 
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We show that the decay mode B — > DrVr is competitive with and complementary to B ri^-r in the search 
for charged-Higgs effects. Updating the relevant form factors, we find that the differential distribution in the 
decay chain B — + Di'-rT~[-^ n~h'r] excellently discriminates between Standard-Model and charged-Higgs 
contributions. By measuring the D and tt~ energies and the angle between the D and tt" three-momenta one 
can determine the effective charged-Higgs coupling including a possible CP-violating phase. 



INTRODUCTION 

The B factories BABAR and BELLE have accumulated 
enough statistics to probe extensions of the Higgs sector of 
the Standard Model. Notably, the decay S+ t+v^ al- 
lows us to place useful constraints on the parameters tan/? 
and Mh+ of the two-Higgs-doublet model (2HDM) of type 
II llll. Here tan [3 is the ratio of the two Higgs vacuum ex- 
pectation values and Mh+ is the mass of the physical charged 
Higgs boson in the model. Since the couplings of iJ+ 
to b's and r's grow with tan/5, 5+ t^Vt probes large 
values of tan (3. Earlier (but less powerful) constraints on 
the 2HDM were obtained by the OPAL collaboration, which 
found tan^/M^+ < 0.53GeV-i from 6(5 XtV^) [2] 
and tan/3/Mjj+ < 0.78GeV-i from B{t iw^,Ur) [3] 
at the 95% CL. The direct search for a charged Higgs bo- 
son through t bH^ at the Tevatron has yielded slightly 
stronger bounds: Mjj+ > 125 GeV for tan /5 = 50 and 
Mh+ > 150 GeV for tan/3 = 70 f^. In the low and in- 
termediate tan/3 regions, the most constraining bound cur- 
rently comes from the FCNC-induced process b 57, which 
yields Mh+ > 295 GeV independently of tan/? (3\. At tree- 
level the Higgs sector of the Minimal Supersymmetric Stan- 
dai-d Model (MSSM) coincides with the type-II 2HDM. The 
coupling of to fermions can be modified by a factor of 
order one due to tan /3-enhanced radiative corrections |6, 7], 
yet this introduces only a few additional supersymmetric pa- 
rameters and the access to the Higgs sector in (semi-)leptonic 
B decays is not obfuscated like in many other modes, such 
as the loop-induced b sj decay. This explains the great 
theoretical interest in the experimental ranges for B{B^ 

T+Ur) ISJ. 

The decay B — > Dtv^ provides an alternative route to 
charged-Higgs effects [9-15]. As we will show in the fol- 
lowing, this mode is not only competitive with 5+ t^Vt, 
but also opens the door to a potential CP-violating phase in 
the Yukawa couplings of the to b and t. B ^ Dtvt 
compares to S ^ rvr as follows: 

i) B{B Dtvt) exceeds B{B tvt) by roughly a factor 
of 50 in the Standard Model. 

ii) B — > Drvr involves the well-known element Vcb of 



the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The un- 
certainty on \ Vub \ entering B rvr is much larger 

iii) B{B TVr) is proportional to two powers of the 
B decay constant fs, which must be obtained with non- 
perturbative methods. Current lattice gauge theory computa- 
tions are struggling with chiral logarithms and fg can only 
be determined with an uncertainty of 30% or more iflill . 
B — > DtVt involves two form factors, one of which can 
be measured in S ^ Dlvg {I ~ e,/i) decays lIlTl llSll . 
The other one is tightly constrained by Heavy Quark Effec- 
tive Theory (HQET) [19-22], so that hadronic uncertainties 
can be reduced to well below 10% once the measurement of 
B Dlvi is improved. 

iv) Unlike B rvr the three-body decay B DrVr 
permits the study of decay distributions which discriminate 
between W+ and H+ exchange jgi fTllflll . 



v) The Standard Model (SM) contribution to S ^ tv^- is 
(mildly) helicity-suppressed, so that the sensitivity of B{B 
TVr) to is enhanced. For B DrVr a similar effect 
only occurs near the kinematic endpoint, where the D moves 
slowly in the B rest frame While the transversely polar- 
ized contribution suffers from a P-wave suppression, the 
virtual recoils against the D meson in an unsuppressed 
S-wave. 

Items iv) and v) strongly suggest to study differential decay 
distributions in B ^ DtVt- The r in the final state poses an 
experimental challenge, because it does not travel far enough 
for a displaced vertex and its decay involves at least one more 
neutrino. In particular, the t polarization, known as a charged- 
Higgs analyzer ifioll . is not directly accessible to experiment. 
To our knowledge, the only theory papers which address the 
uestion of the missing information on the t momentum are 
inil . where a study of the D meson energy spectrum is 
proposed. Another straightforward way to deal with the miss- 
ing information on the r kinematics, which in addition retains 
information on the r polarization, is to consider the full de- 
cay chain down to the detectable particles stemming from the 
T. We have studied the decays t^^Vt, ^ P^Vt, 

and l^viv-r and assessed the sensitivity of the de- 

cay distributions to i7+ effects. We find that the decay chain 
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B DPt-t~[-^ '^~Vt] discriminates between W+ and _ff+ 
exchange in an excellent way. In this Letter we present the 
results for the differential decay rate as a function of the D 
and vr^ energies and the angle between the D and tt^ three- 
momenta for this decay chain. Our result greatly facilitates 
the determination of the effective coupling gs governing iJ+ 
exchange, including a potential complex phase, if e.g. a max- 
imum likelihood fit of the data to the theoretical decay dis- 
tribution given below is employed. A conventional analysis 
combining Monte Carlo simulations of i? ^ Dv^-t^ and 
T~ vr^ i^r decays would be very cumbersome, because the 
B DtVt- differential distributions strongly depend on the 
a-priori unknown value of gs- 



B ^ D FORM FACTORS 

The effective hamiltonian describing B {D)TVr transi- 
tions mediated by or iJ+ reads (with q~u{c)) 



eff 



V2 
mbm 



m. 



- q[gs + gpib] b [t{i - 15)1^^ } + h.c. 



(1) 



The effective coupling constant g p only enters the B rvr 
decay, while B Drvr is only sensitive to gs ■ The _B+ 
meson mass is introduced in Eq. so that B{B rvr) 
vanishes for gp = 1. The above operators as well as mi, are 
defined in the MS scheme. In the MSSM, which is our main 
focus, one has gs = gp- 

The analysis of i? ^ Dnyr requires the knowledge of the 
form factors Fy and Fs which parametrize the vector and 
scalar current matrix elements: 



{D{pD)\crb\B{pB))^Fv{q^ 



1 



Pd 



+ Fs{q')ml^—^qr, 



{D{pD)\cb\B{pB)) = ""1^^ ^^^ Fs{q') , 



nib - nic 



(2) 



where pB and pB denote the meson four-momenta, q = 
Pb — Pd, and r = mpi/mB- It is convenient to intro- 
duce the normalized form factors Vi = Fy 2y/r/ (1 + r) and 
Si = Fs{l + r)/{2y^), as well as the kinematic variable 



w ■ 



(1 + _ q^/ml)/2r. 



(3) 



In the limit of infinitely heavy quark masses itiq — nib, "^c 
(which are properly infrared-subtracted pole masses), both 
Vi (w) and Si (w) reduce to the universal Isgur-Wise function 
^(w), normalized to ^(1) = 1. At the kinematic endpoint 
w = 1, corrections to this limit read 

^l(l) = Vv - ^ {Sr.i + <5i/™q) , 5l(l) = Vv , (4) 

up to 0{a^, l/mg). Here rj^ denotes radiative corrections in 
the limit of equal heavy meson masses, and (5,-ad (^i /mq ) 




FIG. 1: Vector form factor Vi{w). Dots: exp. data fn\ \ with stat. 
errors only. Dashed: fit to parametrization in [22]. Plain: fit to linear 
parametrization in f26il . Dark gray band: form factor with HQET 
constraint at ui = 1; systematic errors dominate at large recoil. Light 
gray band: form factor from HFAG IT^ . 



the first order radiative (l/mg) corrections to the function ^_ 
defined in (t^. The (5i^mQ term depends on the subleading 
function £,3(10 = 1) = Arj{w = 1) and on the HQET param- 
eter A. We take A = 0.5 ± 0.1 GeV, 77(1) = 0.6 ± 0.2 |23ll , 
rjy and Si-^d to 0{as) from Ref. [24], and add a 5% error to the 
form factors at w = 1 to account for higher order corrections. 
We obtain Vi(l) = 1.05 ± 0.08 and S'i(l) = 1.02 ±0.05. 

The semileptonic decay into light leptons B — > D£i'i de- 
pends solely on the vector form factor Vi{w). The mea- 
sured quantity |Vcb 114(11;) was fitted by the BELLE collab- 
oration fll^ to a two-parameter ansatz Vi(w,Vi{l), pi) 1221 
derived from dispersion relations and heavy quark spin sym- 
metry I2M ■ The fitted curve, however, suffers from large sta- 
tistical and systematic uncertainties: |Vcfc|Vi(l) = (4.11 ± 
0.44 ± 0.52)%, pI = 1.12 ± 0.22 ± 0.14 [17]. We thus 
take Vi{l) from HQET instead, use |V'cbJ = (4.17 ± 0.07)% 
from inclusive semileptonic B decays f2y\, and only fix the 
form factor at large recoil w — 1.45 from the data, includ- 
ing the dominant systematic errors in a conservative way: 
\Vcb\ViilA5) = (2.63 ± 0.51)%. The form factor over the 
whole kinematic range is then obtained using a two-parameter 
description Fv{w, , aY), which uses a conformal mapping 
w z(w) resulting in an essentially linear dependence of 
Fv on z 1261] . This linearity in z{w) is confirmed by the fact 
that fitting the B Dlvi data with both Fv parametriza- 
tions without further theoretical constraints essentially gives 
the same result (see Fig. [T]). The sets of parameters corre- 
sponding to the minimal and maximal form factors satisfying 
the HQET constraint at w 1 are displayed in Tab.Ufor both 
parametrizations Vi {w, V\ (1), p^) and Fyiw, Oq , aY). They 
delimit the dark gray area in Fig.[T] We stress that the large er- 
ror band in Fig.[T]at large w is not due to theory uncertainties 
but rather to the large systematic error on |T46|Vi(1.45) from 
il7|]. 

We choose to use only the most recent set of experimen- 
tal data for our numerical analysis. The HFAG |18] treats 
systematic errors in a different way and, including the older 
CLEO and ALEPH data, finds smaller uncertainties at large 
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Parameters 


min. \V^b\F 


max. jVctl-F 


centr. \Vcb\F 


{\V,,\V,{l),pl} 
|V;i,|{af,af}[10-*] 


{0.040, 1.47} 
{0.94, -5.7} 
{1.62,-1.1} 


{0.048, 1.06} 
{1.28, -2.2} 
{2.14, -3.2} 


{0.044, 1.24} 
{1.11,-3.9} 
{1.88,-6.8} 


TABLE I: Parameters {\Vcb\Vi{l), pi} for \Vab\Fv IH and 
{IKfcl^'^, \Vab\aX''^} for \V^i\Fv,s (see |26], = 0, ?? = 2, sub- 
threshold poles: mil-) = 6.337, 6.899, 7.012 GeV and m(0+) = 
6.700, 7.108 GeV lH). Fy is displayed in dark gray in Fig.Q 


Parameters 


min. \Vcb\F 


max. IVctl-F 


centr. \Vcb\F 


{\V.b\V,{l),pl} 
\V.b\{a^ ,aX}[lQ-^] 


{0.038, 1.01} 
{1.03, -1.3} 
{1.78, -5.7} 


{0.047, 1.30} 
{1.17, -4.8} 
{2.00, -7.6} 


{0.042, 1.17} 
{1.10,-3.0} 
{1.89,-6.6} 



TABLE H; Parameters {|V;6|Vi(l), p?} for 

(ii,and{|Ki,|aJ''^, 
light gray in Fig.[T] 



V,b\aX^} for \V,b\F^ 



v,s 



Vcb\Fv from HFAG 
Fv is displayed in 



recoil (see light gray band in Fig.[T] The corresponding mini- 
mal and maximal curves are given in good approximation by 
the parameters in the first two lines of Tab. HI] for w inside 
the B — > Drvr phase space). The vector form factor has 
also beeii studied on the lattice. Computations with quenched 
Wilson i27ll and dynamical staggered ll28ll fermions, however, 
both suffer from potentially large systematic errors, which are 
not fully controlled. In the end, the improvements in the mea- 
surements of the B — > Dlvi and B DtVt modes will 
go together, and | Vcb \ Fy will most likely be best determined 
from experimental data alone. For the time being, we will 
proceed with the conservative estimation of Tab.]!] 

In a similar way, the scalar form factor Fs{w,aQ ,af) is 
constrained by HQET atw = 1, while its value at large recoil 
is fixed from the relation Fs{q^ = 0) = Fv{q^ = 0). The 
resulting parameters are displayed in the third line of Tab.|l](or 
Tab, nil if Fv is taken from MIST ). As expected from the heavy- 
quark limit, the normalized form factor Si is quite close to Vi 
on the whole w range, with slightly smaller errors. 



CHARGED-fflGGS EFFECTS 

The MSSM is a well-motivated new-physics scenario in 
which charged scalar current interactions occur at tree-level. 
Resumming the dominant tan /3-enhanced loop corrections to 
all orders, the couplings gs^p in Eq. ([TJ specify to |13, 30ll 



gs ^ gp = 



(1 + eo tan/3)(l + tan/3) 



(5) 



This particular form holds in MSSM scenarios with Mini- 
mal Flavor Violation (MFV). The loop factor eq arises from 
the quark Yukawa sector and depends on ratios of super- 
particle masses, resulting in a sizable non-decoupling effect 
eotan/3 = 0{1) for tan/5 = C9(50). e,- comprises the cor- 
responding effect for the r lepton. eq and e^. can receive siz- 
able complex phases from the Higgsino mass parameter ^, if 
first-generation sfermions are sufficiently heavy to soften the 



impact of the bounds on electric dipole moments on arg^. 
Beyond MFV also phases from squark mass matrices will eas- 
ily render gs complex. It is therefore mandatory to constrain 
- and eventually measure - both magnitude and phase of gs- 
The type-II 2HDM is recovered by setting — er — 0- 

The B Dtv-t branching ratio has recently been mea- 
sured by the BABAR collaboration llsill : 



^exp 



B{B Drvr) 
B{B Dlivi) 



(41.6 ±n.7± 5.2)% 



(6) 



The normalization to B{B Divg) reduces the dependence 
on the vector form factor Fy and thus tames the main theoret- 
ical uncertainties. In the presence of charged-Higgs contribu- 
tions, the theoretical ratio is approximated to 1 % by 

rl (L544- 



L126 + 0.037 



0.082rs + iVff+) 



10- 0.95 

Nh. =-rcfci?e[gs] (1.038 + 0.076 rs) 
+ r2j5s|2 (0.186 + 0.017 rs). 



(7) 



with TV = {aX /al)/{~ZA), rg = (of /af )/(-3.5), ro = 

(oq /a^)/n, and Vcb = 0.8/(1 — rrTc/mt). The dependence 
V s 

on the slope parameters Ui ' appears to be quite mild. In 
Fig.|2]we compare i?* (right-hand side) as well as B{B 
Tv) (left-hand side) to their one-sigma measurements for 
positive gs and gp. For i?*, we also display the less con- 
servative theoretical prediction obtained from the HFAG vec- 
tor form factor in Tab. HH (light gray band). In particular, we 
obtain the SM estimates 



B{B- 



DrY"^ = (0.71 ± 0.09)% 



and 



S(S° ^ D+T- 



\SM 



= (0.66 + 0.08)%. 



(Error sources: \Vcb\Fv[w), jSj(l), | Kb |). We cannot repro- 
duce the small errors of Ref . lll4|] . 

The B Drvr branching fraction is promising to dis- 
cover - or constrain - charged-Higgs effects, but not to mea- 
sure gs with good precision, as the dependence in Fig. |2] 
is too flat. The differential distribution in the decay chain 
B — > DDtT~ [— *■ tt^i^t] is better suited for that purpose. The 
experimentally accessible quantities are the energies Ep) and 
of the D and tt" mesons, respectively, and the angle 9 be- 
tween the three-momenta po and p^^ . We define these quan- 
tities in the B rest frame, which can be accessed from the 
T{AS) rest frame thanks to full B reconstruction (37]. We 
integrate over the phase space of the two unobserved neutri- 
nos in the final state. Our formulae contain the full spin cor- 
relation between the production and decay of the r, which 
is important to discriminate between SM and charged-Higgs 
contributions. This approach further facilitates the rejection 
of backgrounds from neutral particles escaping detection, as 
in B ^ DD^[-^ 7r^7r°] with an undetected tt'': If the mass 
of the undetected particle is to, this background can be sup- 
pressed by cuts excluding the region around 



cos 6 



(tob ~ Ed - Et, 



2{E 



2[El 



(8) 



FIG. 2: Left: BJB tv) as a function of gp. Light gray band: 6="^ = (1.41 ± 0.43) x lO"* d. Gray band: B*, = (3.86 ± 0.09 ± 
0.47) X 10"^ Isl], main error from B decay constant /s = (216 ± 38)MeV 1 16]. Right: R = B(B DTVr)/B{B Dlvt) as a function 
of gs- Light gray band: 7?="^, see ^ |31]. Dark gray band: i?*, see 0. Gray band: i?"' with the HFAG vector form factor, see Tab.|ll] SM: 
^?*(SS = 0) = 0.31 (dark gray) [ ±0.02 (gray)]. 



We obtain the differential distribution 

dT{I3 Di>-rT-[-f TT-V-r]) 



= G%fi\VudV\VcbVTr 



(9) 



dEo dEjr dcos6 
X [CwiFv.Fs) - Cwh{Fv,Fs) Re[gs] + CHiFs)\gs\''] 

with form-factor-dependent functions of Ed, E^^, and 0086* 
for the SM (Cw), interference (Cwh), and Higgs (Ch) con- 
tributions, given as follows for vanishing ttIt^ (this approxi- 
mation, which is good to 1%, is not used in our numerical 
analysis). 



Cw = K 



2 -l 
P{P-p.f 2{P-l){P-p^) 



Cwh = 2k, 



Ch — n 1 



ip-rr ■ ly Ptt ■ I 

4 (1 - r')Fs 



1 - mc/mb 

5 (1- 



(1 - nic/mby 



P-l- 



1 



(36-1) 



(10) 



2ptt ■ I 



where nic and rub must be evaluated at the same scale so that 
mc/mb = 0.20 ± 0.02 133(], and 

P = FvipB +Pd) - {Fv - Fs)^^. -{pb~Pd) , 



I 



128 7r4 iTiBmr 
Pb - Pd ~ Ptt ■ 



Ptt ■ I 



2p^ ■ I 



q 



{pb -Pd) 



(11) 



The dot products appearing in Eqs. ( fTOl i and ( fTTT i are re- 
lated to the energies, momenta, and the angle 9 measured 
in the B rest frame as pB ■ I ~ mB{mB — Ed — E-^), 



PD-l ^ EnimB - Ed - E^) + \pi 



IpdIEttCos 0, 



Pn-l = Ej,{mB ~ ED) + \pD\Ejr COS 9, and pbPd = msED- 
Further = (290.6 ± 1.0) x 10~^^s is the r lepton Hfetime, 
/tt (130.7 ± 0.1 ± 0.36) MeV the pion decay constant, and 
the CKM matrix elements are = 0.97377± 0.00027and 
\Vcb\ = (41.7 ± 0.7) X 10"^, the latter being well determined 
from inclusive semileptonic B decays i25ll . Remarkably, one 
can probe a CP-violating phase of gg by exploiting the shape 



of the distribution in Eq. (|9]), which is not possible from the 
branching fraction of either B DtVt^ or B ^ tv^-- 

For illustration, we show the differential decay distribution 
including charged-Higgs effects in comparison with the SM 
for the meson energies Ed = 2 GeV and E^ = 1 GeV, so that 
the whole range of cos 9 is kinematically accessible. In this 
particular region of phase space the SM rate is strongly sup- 
pressed for cos6' = —1. For a large scalar coupling gs — 2 
(Fig- El left), the Higgs contribution dominates the rate at this 
point (dark gray band), so that we can clearly distinguish it 
from the SM (light gray band). The experimental informa- 
tion from B{B TVr) constrains |1 — gp\. For real gp 
this permits a range near gp — Q and another range around 
gp — 2. In the MSSM situation with gp = gs, the case 
gs — 2 therefore is in agreement with B tVt, but can be 
confirmed or ruled out by measuring our distribution. The dis- 
crimination potential for the phase of gs shows up in the light 
gray band: It coiTesponds to a complex gs = 1 -\- i, which 
yields the same B{B tVt) as gs — 0, 2. The B Drvr 
branching ratio alone may also help to distinguish between 
these solutions, depending on the future experimental value 
of B{B DtVt), see Fig.|2l For general gs values, a fit to 
the triple differential distribution in Eq. (|9]l would excellently 
quantify charged-Higgs effects, especially once better exper- 
imental information on the form factors is available, as we 
illustrate with fixed D and vr^ energies in Fig. [3] (right-hand 
side) for gs — 0.5. Such a fit would combine information 
from different parts of the phase space, and thus resolve much 
smaller gs values. A more precise quantitative analysis would 
require the fit to actual data, and thus goes beyond the scope 
of this paper. Still, keep in mind that even with more precise 
B — > tVt experimental data and improved estimates of Jb 
and \ Vub\, a value of(7p~0.2 — 0.3 will be very difficult to 
exclude with B{B tv^)- B DrVr is thus definitely 
competitive. 

As mentioned in the Introduction, a similar analysis was 
performed for the other r decay channels ^ p^Vr and 
t^viVr, which together with ir^Vr constitute 

more than 70% of the r branching fraction. Ultimately, a com- 
bined analysis of all these modes is desirable in order to ex- 
ploit the available and forthcoming experimental data in an 
optimal way. 
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FIG. 3: -> D+UrT~[-^ Tv'Ur] angular distribution for Ed = 2GeV and = IGeV. Left: gs = 0, 1 + i,2. Riglit: gs = 0,0.5 
(dark gray: witliout uncertainties in Fv (w) and Vet, errors from 5*1 (1) and rfic/rnt). Tlie conservative form factor estimates of Tab. Uwere 
considered. 



CONCLUSIONS 

We have studied charged-Higgs effects in a differential dis- 
tribution of the decay chain B — > Di>TT~[^ tt^Vt], which 
has the following advantages over the branching fractions 
B{B TVr) and B{B Drvr): 

i) The Higgs coupling constant gs can be determined from 
the shape of the distribution in sensitive phase space regions. 
This analysis should be possible with current B factory data. 

ii) The dependence on both \gs\ and Re [gs] allows to quan- 
tify a possible CP-violating phase. Since our decay distribu- 
tion is a CP-conserving quantity, the phase of gs is deter- 
mined with a two-fold ambiguity. In the MSSM such a phase 
stems from the /i parameter or the soft breaking terms and 
enters through tan /3-enhanced loop factors. B Dtv com- 
plements collider studies of these phases ll34ll . 

The main uncertainties stem from the form factors. One 
can gain a much better accuracy with better data on the vec- 
tor form factor Fy- The recent B Dlv^ measurement by 
BABAR 1 35] furnishes promising data for a new fit. 

Within the MSSM, one will be able to place new constraints 
on the tan (3 — Mh+ plane, once our results are confronted 
with actual data from the B factories. If tan/3/M^+ is in- 
deed large, there is a fair chance to reveal charged-Higgs ef- 
fects ahead of the LHC. 
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